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Abstract 

A formal linear vector field representation for scientific equations is developed 
to rationalize the intuitive methods that are constantly employed. It is shown 
that unlike mechanical units that appear in the basis of the space, the reduced 
temperature and Boltzmann parameter cannot be described by the basis set 
individually and can only be described as a product. Further, the definition 
and determination of temperature is dependent on theory and not on standard 
mechanical units. It is shown that there is no reason to reduce the number 
of degrees of freedom in temperature determination via equipartition since 
stochastic variables are involved, and this observation is significant in that the 
temperature variable reported in simulation studies would have a discrepancy 
to the extent of using the decreased number of freedom, which is most cases 
is not large nor significant. The standard assignments used in reduced units 
do not lead to errors because operationally the resulting reduced temperature 
parameter represents the reduced product of the temperature and Boltzmann 
parameters. The non-independence of these quantities explains why entropy 
and other associated functions cannot be calculated directly, but are always 
scaled in dimensionless increments of the Boltzmann parameter 

1 Introduction to vector field properties of 
scientific equations 

A scientific result of magnitude Cj may be represented as Cj. [11™ i M" 1 ] where 
usually a.i G Z, q G F for field F and where normally the restriction a € 
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K is utilized only but there is no reason why q or cti cannot belong to a 
larger (e.g. complex) field. The Mj are the so-called fundamental units; 
e.g. c. [M a L^T 7 ] may describe a quantity of magnitude c characterized by 
exponents (a,/3,7) in mass (M), length (L) and time (T) respectively. The 
set {Mj}, in general, and in particular M, L, T is isomorphous to any fixed 
members in F where the exponents of products (operator ©) is concerned, i.e. 

[UT= i M i Ql ] ® [UT=1 M i ft ] = [UT=i Mi^ +ft )] . Let d a = []JT=i Mi a ' ] define 
the " 'unit dimension" ' or basis of the scientific quantity; then d a is uniquely 
characterized by the m-tuple (a±, 02, ... a m ). Define the zero element and 
unit element of d Q as 0. [U7=i M i a ] and MldXiW* 1 ] respectively. We 
define the addition operator © such that 

The scalar multiplication of element c ai .d a by q € F is defined as 

q(c a .d a ) = (c ai .q).d a . (2) 

For any vectors a, /?,7 where a\ = aj.d a in space d a , the operator is 
commutative, associative, and properties (1) and (2) ensures that = 0.d Q 
and Qai = — cti.d a so that a © = a for all a in d a and a © (0a) = 0. 
From Q and (J2J), we infer that {a} forms a one dimensional vector space. 
From experience, only similar quantities may be added, which leads to Axiom 
1. 

Axiom 1 In scientific equations, the addition operation between vectors be- 
longing to different dimensional basis is not possible, so that if d a ^ dp then 
oii.d a © Pi.dp is not defined. 

From Axiom ^ it follows that c a .d a is a 1-Dimensional vector space with 
the unique zero 0.[J|™ 1 Mj Qi ] which is not equatable with the zero vector of 
another space cg.d^ . In this sense, the dimensional bases d 7i are orthogonal 
to each other. 

Definition 1 The product operator © is a mapping P such that 

P : {c a .d a , cp.dp} -> Ca.cp. [da x d^] (3) 

and is defined as follows: [c a .d a } A ®[c p .dp\ B = (c a ) i (c/[nr = iMi (Aa,+Rft) 

where [d a x d^] = n^Li Mi (AQl+R/3l) , and A,B G F (but commonly re- 
stricted to 7j).The properties 0/© operator are as follows : 

1. if B above is negative, then cp 7^ 0. 

2. it is symmetric, i.e. [c a .d a ] A © [cg.d j g]' B = [c^.d^]" 8 © [c Q .d a ] A . 

3. it is associative, where for all elements II, 0,£, the following obtains 

(n©e)©£ = n©(e<g>£). 
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4- it is distributive over addition , i.e. 



[c a .d a ] A [cfo.dp © Cfo.dp] = Ca-Cfa. [d a A x d^] c A .cp 2 . [d a A x dp] 

= ^.c ft + ^.c ft ).[d/xd^]. (4) 

5. the divisor operator is defined such that 0c a .d a = 0(c a ) _1 [ YYiLi M[~ ai ] . 

Definition 2 A general mapping transformation i^ Q '^'■■■ 7 which maps the 
domain (c a .d a cp.dp, . . . c 7 .d 7 ) to a range in d^, CR.dR in a series of opera- 
tions involving only the © and operators and scalar multiples is defined to 
be consistent. 



Axiom 2 Scientific equations are consistent mappings. 

From Definitions (EK^; consistent mappings are of the form 



a, (3,... 7 



(c a .d a , cp.dp . . . c 7 .d 7 ) = fi(c a , op, ... c 7 ). 



(a, A -7) 



i=l 



(5) 



where each a, (a = a, /3, . . . 7) are m-tuples a = (o~i, 02, 



Remark 1 (a) The mapping here involves the mapping of a collection of 
1-D vector spaces {Wi} — > V where the domain V is also 1-D with basis 
d v where d v = [\\T=i M ^]- 

(b) Transformation |2J) defines a scientific equation where fi(c a , cp, . . . c 7 ) E 
F a member of the scalar field, representing a physical quantity with di- 
mensions characterized by the functions f[. 

From (J3J) and ©, the field properties under (.,+) for a scalar function 
fi(a, (3, • • • 7) — ► F is isomorphous to (0, ©) under the 1 — 1 mapping a — > 
a. [M Q ] and so one can write 



/,(«, (3, ...7) [M F '[ f i( a ^-^l =f l (a.[M a ], ...7.M 7 ) 



(6) 



where the (., +) operators of the scalar function / are replaced by the (0, ©) 
operators respectively in the R.H.S. of ©• The above is the reason why some 
people speak of reduced units as being '"unitless"' ^ p. 199]. An example of 
a scientific equation for /, is the simple Lennard-Jones 12-6 potential, 



v LJ (r,e,a) = Aire 



12 



<7\6 

r 



(7) 



where the basis unit vectors are 1. [M a L^T 7 ] where M,LandT are the mass, 
length and time base unit symbols (e.g. kilogram , metre and seconds in S.I. 
units); v LJ is in this notation entirely unitless as all variables are members of 
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the R field. From ©, the transformation function .F' 
the form 

i^f'- 7 (e. [M°L 2 T- 2 ] , r. [M !^ ] . . . a. [M°L 1 T°] 



e,a). [M 1 L 2 T~ 21 



is in this case of 



(8) 



where the units of the potential are characterized by [M 1 L 2 T -2 ]. Scientific 
equations are described relative to its dimensions; for two basis systems M 
and M* describing the same physical phenomena, each having the same span 
length m for its basis vectors , where for system M, the basis is written 
M a M% . . . M7n and similarly for the M* basis, each dimensional symbol is 
related linearly, i.e.M* = AiM;. This leads to Axiom 3. 

Axiom 3 Scientific equations describing phenomena relative to the basis 
l-[n2=i Mi ] may be described relative to another basis l-[II2=i Mi* ] through 
fixed scaling parameters Aj where 



I I A? 



\i=l 



IK 



i=i 



i. 



f[ Mi* a ' 



i=l 



(9) 



The vector spaces are therefore linearly dependent. Scientifically, the choice 
of units cannot led to different physical phenomena through the scientific 
equations describing the system trajectory, which is expressed in Axiom 4. 

Axiom 4 The mappings describing scientific laws in scientific equations are 
independent of chosen unit dimensions, so that for any scientific equation 



P a,/3,...7 



(c a .d a , cp.dp ...Cy.d^ 



pa, /?,... 7 



(c*.d a *, c* p .dp* ... C ;.d 7 *) (10) 



Axiom and © yields for the R.H.S. of (|10[) the following 



a,A-7 („* 



c*.d a *, cp.dp* ...c*.d 7 *) = /j(c*, c*j, ...c*). 



i=l 



Here we write /* = /j. Axiom El and (f X If) implies 



(a;, ft, ...7i) 



Ff'^'"' 7 (c* A Q .d Q , c^A^.d^ ...c*.A 7 d 7 ) 



ft,-7i) 



J^J M. ft,- 7i) 



(12) 



But 

ft... 7 



Ff'^'" 7 (c Q .d a , c^.d^ . . .c r d 7 ) = /j(c a , c/3, . . .c 7 ). 



ft,- 70 



i=l 



4 



(13) 



Thus we have 

r:,...r:;.nAf"" j ' = /(a,, A, • • • 7i)- 
Axiom |21 and © give 



or, 
c* = 



wrr=iM Q / 



A, 



i=l 



i=i 



(14) 



(15) 



(16) 



From the above, the star operator (*)corresponding to a change of unit basis 
l.[M] to l.[M* ] and any scalar function with variables c a linked to the vector 
c a .[M„] may be written 

/ (c a , Cp, . . . Cy) = /(c* , . . . c*) 

C(3 C 7 



/ ( \ • \ • 
with (A Q , A/3 . . . A 7 ) given in l|16jl. 



A, 



-) 



(17) 



2 Discussion and verifications 
2.1 (a)some standard applications 

In the laboratory [M] basis, the interparticle potential has sometimes been 
modeled according to j7}; a change in the unit basis implies converting the 
laboratory frame of units given in (J7J) to another [M* ] with form given by (jllj) . 
In this situation, let there be a unit basis of length such that a. [M^ ] = l.[Mr ] 
or in detail a[M L ] = a. [M°L 1 T°] = l.[M* L ] = 1. [M^L^T* ]. Further, let 
the energy scale as e.[Ms] = l.[M^]. In @, c and r are linked to [Ml]; the 
transforming operator on v LJ * = v LJ (a*,e*,r*) yields 
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V LJ * = 4vre 




(18) 

which is a reduced potential used in simulations. In the above case, the 
bases for energy and mass are [M E ] = [M 1 L 2 T~ 2 ] and [M M ] = [M^T ] 
respectively. If 3 scales are chosen for the A's, e.g. (e, a, m), then other 
quantities are fixed relative to it. From the kinetic energy, we allow e = ^mv 2 , 
then the scaling for the velocity is also determined, and so v 2 is a fix quantity. 
We require v' , the A parameter for velocity scaling. Applying Q17|). the results 
are e* = 1 = \(^r) 2 or v' 2 = \ = ~ which results in the scaled velocity 
v* = in accordance with standard results 0. 
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2.2 (b)the temperature parameter 

The above theory represents scientific numbers in terms of c a . [M ai L a2 T a3 ], 
and clearly, if a unit cannot be expressed in terms of the defined basis given 
here, then absolutely no scaling parameter A exists for that quantity, and 
in particular (a.c)* ^ a* .c* since the isomorphic properties of the above- 
mentioned operators only apply to quantities with an associated dimension 
c .[M Q ]; for example scaling with fixed parameters (m,o~,e) means that the 
velocity scaling parameter A^ is defined such that v* = j- is (^) 1//2 , where 
the dimension of this quantity, as well as the scaling parameters all have the 
form c a .[M Q ]. Temperature (and other associated properties) do not emerge 
directly from mechanical analysis - where the basis dimensions are defined - 
but from theoretical models or definitions, and in particular cannot be defined 
in the form c Q .[M Q ] . It is therefore misleading to suppose that a Xt divisor 
must exist for this quantity such that T* = which is often suggested e.g. [2J 
Table 1.]; the actual implication, as shown below is the definition of another 

temperature scale T = ksT where T* = y^r~)- Likewise the Boltzmann 
coupling parameter ks has a unit which is reciprocal to temperature, imply- 
ing that it cannot be expressed in the form cg.fM^ ] and hence this parameter 
cannot be reduced in isolation by setting "'/cb = 1> so that the MD unit of 
temperature is also defined"' [HI pp. 15-16]. This assumption of the funda- 
mental mechanical autonomy of the coupling parameter and temperature is 
however rather standard and pervasive p. 200] and in nearly all cases, the 
standard assignments are correct, which is explained below. 

In MD simulations, the temperature is determined via the classical equipar- 
tition theorem (which is known from experience and quantum mechanics not 
to obtain at lower temperatures for the free vibrational and rotational modes) 
for all temperatures from the mean kinetic energy of translation where 



and where D = 2, 3 for two and three dimensional systems respectively for 
an N particle system where the angle brackets denotes some chosen aver- 
aging algorithm that is thought to approximate the outcome if the P den- 
sity function mentioned below is used in an exact evaluation. Investigators 
claim that due to "'conservation of momentum'", there are D degrees of 
freedom that must not be counted in ()19f) .[5l HI p. 16 and pp. 46-47 respec- 
tively]. The probability density function P over the (p, q) momentum and 
positional coordinates for the above averaging process is P = exp / Q 

with Q = J dv exp ~^^ q ' ) dp dq defined as the phase integral, which is the 
analog of the quantum partition function for the canonical distribution. W is 
the Hamiltonian , written in most classical simulations with the momentum 
and potential coordinates separated viz. H = + Using this P 




(19) 
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density, each particle yields < pf/2m >= y(&gT) , where clearly the defini- 
tion of T is the consequence of stochastic averaging; it is not a fundamental 
mechanical quantity associated with a basis dimension. Extending this result 
to N particles is thought to yield ()19|) , which is routinely used to compute 
the temperature in simulations. Averaging the kinetic energy of TC using P 
(which already has the constraint over the entire ensemble of total energy 
conservation - from where the (3 parameter from the Lagrange multiplier rep- 
resents the ksT term and the a term from the Lagrange multiplier refers to 
total particle conservation which is featured in the phase integral or partition 
function) does not indicate any such reduction of degrees of freedom. Indeed, 
from the Gibbs' postulate of the equivalence of ensemble average to that of the 
time average of a particular system, one can view each particle as a system, 
from which we can expect that the time average of the mean kinetic energy 
of any given particle would equal ^ksT with the same ksT as given in (|19|) 
for the entire system of which it is a part. Further, the supposed reduction of 
the number of degrees of freedom implies that no temperature can exist for 
a system comprising of D or less particles, and this is patently absurd, for it 
is eminently feasible to conceive and implement thermostats in MD for such 
systems. Hence the proper form of the temperature must be derived from the 
probability density function and/or Gibbs' postulate with the average energy 
per particle given as S&bT to yield 



M = ^ W (20) 




if the classical kinetic energy is used as an indicator. The error reported in 
studies would be due to the degree of difference between (|19j) and (|2()j). 

From the vector space calculus here, (ksT) can be expressed as e^.pVLj] 
, but (ksT)* ^ k* B T* since the temperature and coupling constant is not 
separately definable. Hence another temperature parameter T may be de- 
fined where T = ksT, so that (T)* = (ksT)* and hence T* = (ksT)* = 

Cwd (Sill m *X* 2 )) an d in particular T* = 0^r-\ which is the standard 
assignment [2], which leads to standard and consistent results, provided it 
is understood that T = k^T and not T = T, which is the normal under- 
standing, where it is assumed [5] that k* B = 1; even with this unfortunate 
assumption, together with the autonomy of variables, the correct results are 
derived because of the following equation sequence: k* B T* = (ksT)* = or 
T* = ^jpr- = which is the fortuitously correct result with k B = 1 but with 
the incorrect algebraic assumptions, since kg is not independent, it cannot be 
arbitrarily set to a value. 

There are clear-cut consequences that follow from whether (a) an indepen- 
dently scaled k* B exists, or, (b) where this is not the case. In the simulation of 
entropic quantities based on the Boltzmann postulate for entropy S given by 
5 = ks In W , case (a) implies S* = \a.W, or a direct determination is possi- 
ble for the entropy; (b) suggests a work-around, such as scaling equations as 
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(i)5" = (S/ks) or (ii) E 1 = (TS) in constant temperature studies. The tech- 
niques used in entropic studies use(i) and (ii) or a variant of (ii) through the 
determination first of the chemical potential scaled in /cbT units ^ pp. 246- 
249], or in coupling methods (utilizing thermodynamical integration), S/ks is 
the variable that is scaled in simulations (method (i) ) 1, p. 260]. It is proposed 
here that simulations utilize variants of either (i) or (ii) above as a direct con- 
sequence of (b), that an independent ks and T (or k* B and T*) does not exist 
for scaling. 
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